Hydrodynamic behavior in expanding thermal clouds of ®^Rb 
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We study hydrodynamic behavior in expanding thermal clouds of *^Rb released from an elongated 
trap. At our highest densities the mean free path is smaller than the radial size of the cloud. After 
release the clouds expand anisotropically. The cloud temperature drops by as much as 30 %. This is 
attributed to isentropic cooling during the early stages of the expansion. We present an analytical 
model to describe the expansion and to estimate the cooling. Important consequences for time-of- 
flight thermometry are discussed. 

PACS numbers: 03.75.Hh, 03.75.Kk 



I. INTRODUCTION 

The anisotropic expansion of a condensate after re- 
lease from a trap is one of the best known features of 
the Bose- Einstein condensed state 0, Q • The anisotropy 
arises because the condensate expands most rapidly in 
directions where it was originally most confined. The 
interest in this phenomenon is further growing, in par- 
ticular since the observation of anisotropic expansions in 
non-condensed Bose gases 0, 0| and in degenerate Fermi 
gases 

Anisotropic expansions are indicative for hydrody- 
namic behavior. It is well known that Thomas-Fermi con- 
densates can be described by the classical Euler equation 
for potential flow of a non- viscous gas Q. Therefore, they 
behave hydrodynamically even at very low densities. For 
classical clouds the situation is density dependent. At low 
densities, where the mean free path is large compared to 
the size of the cloud (coUisionless regime) , the expansion 
proceeds under free flow conditions (free expansion). The 
motion of the individual atoms is described by a single- 
particle Hamiltonian and the expansion is isotropic. Re- 
ducing the mean free path to a value smaller than the 
dimension of the cloud allows the introduction of a hy- 
drodynamic flcld and leads to a crossover to hydrody- 
namic behavior (hydrodynamic expansion). Little differ- 
ence is to be expected between the expansion of a con- 
densate and that of a fully hydrodynamic thermal cloud 
0. Both the coUisionless and the hydrod yna mic regime 
were studied theoretically (see H, 0, 0, ^| and refer- 
ences therein). Also the influence of mean-fleld effects 
[T^ and the crossover between the two regimes were an- 
alyzed theoretically |0| and numerically tIJ. 

It is important to understand the crossover to hydrody- 
namic behavior in thermal clouds. From the fundamen- 
tal point of view it is important to quantify the hydro- 
dynamic properties as these affect the coupling between 
condensates and thermal clouds. From the experimental 
point of view it is vital for the correct interpretation of 
timc-of- flight absorption images of dense atomic clouds. 
Previously the crossover regime in thermal clouds was 
probed in experiments at MIT with a dense gas of ^^Na 



atoms ^3 ^-iid at ENS using cold metastable triplet "'He 
|l6j . In Amsterdam the crossover regime was observed 
in experiments with ^^Rb Very pronounced hydro- 
dynamic conditions were recently reached by exploiting a 
Feschbach resonance in fermionic gases 0, 0, 0, 0, ^3 ■ 
Hydrodynamic behavior as observed in collective excita- 
tions is reviewed in refs. 0,113 ■ 

In this paper we focus on hydrodynamic behavior as 
observed in the expansion of dense thermal clouds of 
^^Rb, extending a brief analysis presented earlier in the 
context of the BEG formation experiments in Amster- 
dam U . The clouds are prepared in an elongated trap at 
a temperature Tq, just above the critical temperature for 
Bose-Einstein condensation. At the highest densities the 
mean free path is less than the radial size of the cloud. 
After release from the trap the clouds expand anisotrop- 
ically and their temperature drops by as much as 30%. 
The behavior is intermediate between that expected for 
coUisionless clouds, where cooling is absent, and pure hy- 
drodynamic behavior, where the gas cools to vanishing 
temperatures. 

We show that the expansion in axial direction is similar 
to that of a coUisionless cloud at a temperature < Tq. 
This 'axial' temperature can be identifled with the tem- 
perature reached at the moment when the expansion 
ceases to be hydrodynamic and the cooling stops. Ra- 
dially, the expansion proceeds faster than that expected 
for a coUisionless cloud and can be characterized by a 'ra- 
dial' temperature Tp > Tq. For our conditions, the mean 
fleld of elastic interaction contributes ~ 20% to the total 
energy in the trap center. We show that this only has a 
minor effect (3%) on the expansion behavior. The con- 
sequences for time-of-flight thermometry are discussed. 



II. EXPERIMENT 

In our experiments we load a magneto-optical trap 
with approximately lO'^'^ atoms from the source described 
in 21]. After optical pumping to the \Si/2, F = 2,mF = 
2) state typically 4 x 10^ atoms are captured in a loffe- 
Pritchard quadrupole magnetic trap. Then the gas is 
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compressed and evaporatively cooled to a temperature 
just above Tc- The radio-frequency (rf) evaporation is 
forced at a final rate oi v — —433 kHz/s down to a value 
vi = 740 kHz, that is 120 kHz above the trap mini- 
mum Bq = 88.6(1) /iT as calibrated using atom laser 
output coupling 1221 • As the final ramp down rate is 
— z>/(t'j^ — i/p) ~ 4 s~^, i.e. slow compared to both ax- 
ial and radial trap frequencies uu^ — 2tt x 20.8(1) s^^ 
and Wp = 2-77 X 477(2) s~^, the evaporation proceeds 
quasi-statically and yields a sample characterized by a 
sing le uniform temperature and an equilibrium shape 
|23j. The preparation procedure is completed by 20 
ms of plain evaporation at rf- frequency vi. This pro- 
cedure leaves us with N = 3.5(3) x 10^ atoms at density 
no = 3.6(6) X lO^'' cm~'^ in the trap center and temper- 
ature To = 1.17(5) /iK. 

A. Knudsen criterion 

To establish the collisional regime we calculate the 
mean free path and the atomic collision rate. The mean 
free path in the trap center is given by the usual expres- 
sion for a uniform gas |25j | at density ng, 

Ao = — — w3pim, (1) 
v2rioCT 

where a — 87ra^ is the clastic scattering cross-section in 
the s-wave limit with a — 98.98(4)ao the scattering length 
[2^ . The atomic collision rate in the trap center is |25ij 

= \/2noUt,,(7 » 6000 s -\ (2) 

with vth = {SkBTo/nmY^'^ as the thermal velocity. 

The gas behaves as a hydrodynamic fluid if the mean 
free path is much smaller than the relevant sample size 
(Knudsen criterion). Defining the axial (l^) and the ra- 
dial (Ip) size parameters of the density profile in a har- 
monic trap, see Eq. , the Knudsen criterion can be ex- 
pressed as 

^ ^ C.,Te « 1, (3) 
H 

with i g {p, z}. For the axial direction the Knudsen cri- 
terion is very well satisfied, llIzTc ~ 0.02. For the radial 
direction we calculate ujpTc « 0.5. In this direction we 
operate in the middle of the crossover range between the 
collisionless and hydrodynamic regimes. 

B. Time-of-flight analysis 

In the crossover between hydrodynamic and collision- 
less conditions the time-of-flight analysis is non-trivial. 
Unlike in fully collisionless clouds, the velocity of the in- 
dividual atoms is not conserved because the gas cools as 
it expands. Unlike in fully hydrodynamic clouds, cool- 
ing will only proceed during a finite period. Obviously, 
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FIG. 1: Expansion measurements for (a) axial and (b) ra- 
dial direction. The error bars represent two standard de- 
viations. The solid lines represent Eqs. Q and with 
1^{U) = 116 ^m, = 0.83 /iK and Tp = 1.35 /iK. Note the 
difference in vertical scale for the two panels. The dashed line 
represents the asymptotic expansion behavior in axial direc- 
tion. As the initial radial size is very small the radial expan- 
sion is already asymptotic by the time the first data point is 
taken. 

if the temperature drops during the expansion the ques- 
tion arises how to properly extract the temperature of 
the cloud from a time-of-flight absorption measurement. 

In Fig. n we plot the measured axial and radial cloud 
sizes, Izit) and lp{t), as a function of expansion time t. All 
data were collected during a single run within 2.5 hours, 
keeping track of some drift in the ofset field (22l. Each 
data point corresponds to the average of about 20 mea- 
surements, with the error bars representing the standard 
deviation, typically 2% of the average value. The cloud 
sizes were determined with the usual procedure (see for 
instance ref. (l^), i.e., the expression for the column den- 
sity of an ideal Bose gas trapped in a harmonic potential 

n2{z,p) = n2og2[De-^'-^''^'^^"-^P^'^^'^^']/g,[D] (4) 

is fitted, after transformation to optical density, to the 
images [22]. With this procedure we obtain values for 
the sizes lz{t) and lp{t), the degeneracy parameter (fu- 
gacity) D and the peak column density 7120 |33|. We use 
the notation ga[x] = X^z^i^V'"- The fugacity provides 
together with the initial sizes a self-calibrating method 
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for the total atom number provided the average trap fre- 
quency w = {ujpLJzY/'^ is known, 

N^9m{^ymitio)- (5) 

In practice only the axial size lz{Q) is used because the 
aspect ratio is accurately known. The measured peak 
column density n2o is not used in our analysis [3ll |. 

Due to the presence of the elastic interactions between 
the atoms the density distribution will be slightly broad- 
ened and deformed l33 | . Calculating the variance of 
the distribution (z^) using the recursive expression for 
the density to first order in mean field J7mf(r) — 2gn{r), 
leads to 

^muj^jl(Q) ~ kTo + E,-ni, (6) 

where £',„f = g f n^(r)dr/ J n(r)dr is the trap averaged 
interaction energy with g — {AiTh^/m)a the interaction 
coupling constant .33]. The variance was related to 
the size parameter using (z^) = iZ^(0)g4[£']/(ji3[I?] [33| . 
Equivalently, treating the mean field as an effective po- 
tential we may write 

\mCjlll{Q) = fcTo, (7) 

where ujz represents a 'dressed' trap frequency that re- 
produces, for an ideal gas at temperature Tq, the same 
cloud size, 

ul=ul{l-il (8) 

where ^ = E,,,f/{kTo + E^i) w 0.03 |3|. 

To describe the expansion behavior analytically we in- 
troduce a schematic model in which the expansion is 
treated as purely hydrodynamic up to time t — and as 
purely coUisionless beyond this point. At the density 
has dropped to the level that no further collisions take 
place and the atomic velocities remain frozen. The axial 
expansion is represented by 

Iz (t) ~ [ll iU) + {2kBTjm) {t - U) 2]i/2. (9) 

The presence of t^, slightly shifts the asymptote of the 
expansion curve. The radial expansion is asymptotic for 
all times relevant in the experiment, 

lpit)^[2kBTp/m]'/H. (10) 

In this case the shift of the asymptote is negligible. The 
parameters and Tp represent apparent axial and radial 
temperatures corresponding to the asymptotic expansion 
velocities of the cloud in both directions, 

= lim (t) = (2fcsr,/m)i/2, (11) 

t — >oo 

with i g {p,z}. Note that Eqs.® and ((101) reduce to 
the usual expressions for isotropic expansion of fully col- 
lisionless thermal clouds in the absence of a mean field 
when with Tz=Tp = Tq (see e.g. ref.jH). 



III. RESULTS 

Fitting Eq.Q to our data, the degeneracy parame- 
ter was verified to be constant during the expansion to 
within experimental error, D — 0.95(4). Once this was 
established we determined the cloud sizes by refitting 
all data with a fixed value D = 0.95. The results are 
shown in Fig.^ (solid bars). Fitting Eq.@ to the re- 
sults for the axial sizes we obtain the initial axial size 
lz{0) — lz{t*) = 116(2) ^m and the 'axial temperature' 
Tz — 0.83(4) /iK. The fit is shown as the solid line in 
Fig-ffli aiid is insensitive to any reasonable choice of i*. 
Fitting Ea. (|10|) to the radial data we obtain the solid 
line in Fig.^, which corresponds to Tp = 1.35(6) /iK 
[ssj . For all these results statistical errors are negligi- 
ble. The quoted errors represent the uncertainty in the 
determination of the fugacity. 

From the initial axial size we calculate with Eq.jTI) 
To = 1.17(5) fiK. Then, the central density no = 3.6(6) x 
10^'* cm"-^ follows with 

no = g3/2[D]/Al (12) 

where Ao = [27111^ /mkTo]^^^ is the thermal wavelength 
at temperature Tq. Using Eq.® to account for the 
mean field broadening we calculate with Eq.© N = 
3.5(3) X 10^ atoms. The error bar reflects the strict condi- 
tions on the atom number imposed by a known fugacity. 
We return to systematic errors in the section on ther- 
mometry. 

The results presented here indicate a slightly decel- 
erated expansion in axial direction, Tz/Tq — 0.71(2), 
and a slightly accelerated expansion in radial direction, 
Tp/To = 1.15(3). This corresponds to an 'inversion' of 
the aspect ratio, which is demonstrated in Fig.|2]by plot- 
ting the aspect ratios for the same data set as used in 
Fig-ffl At t = 12 ms of expansion the cloud shape crosses 
over from a cigar-shape to a pancake-shape. The solid 
line represents a fit to the expansion model to be dis- 
cussed below. 

For coUisionless samples the expansion is expected to 
be isotropic. This was verified by reducing the density 
by a factor of 30 (open circles in Fig.|2Jl. In this case 
the expansion is indeed isotropic (dashed line), Sz/sp = 
1.02(4). 



IV. EXPANSION MODEL 

To interpret our results for To, Tp and Tz we divide the 
expansion in two stages. During the first stage (i < t*) 
the expansion is treated as purely hydrodynamic and is 
described by scaling theory |3, 0| . All data are taken 
during the second stage {t > t^,) for which the expansion 
is treated as coUisionless. 
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FIG. 2: Aspect ratio of a hydrodynamically expanding cloud 
as function of expansion time. The error bars represent two 
standard deviations. The change from a cigar- to a pancake- 
like shape is evident as the data points cross the value of 
Ip/lz ~ 1. The open circles represent low-density clouds ex- 
panding isotropically. The solid and dashed lines represent 
fits of Eq. 1)35^ to the data. 



A. Hydrodynamic stage 

During the hydrodynamic stage (i < t,) wc treat the 
expansion as isentropic, i.e., the gas cools while convert- 
ing random motion into directed motion just as in the 
supersonic expansion of an atomic beam [s^ . As for 
isentro pic expansions the degeneracy parameter D is con- 
served we find, using Ea. fT^ . that the temperature 
decreases according to 



(13) 



Turning to scaled size parameters, bi {t) = li{t)/li{Q) 
with i g {p, z}, the density ratio is conveniently written 
as 



1 



no bl{t)h,{ty 



(14) 



We note that for our elongated clouds (w^/wp <C 1) the 
axial size remains practically unchanged during the early 
stages of the expansion. Therefore, setting 6^ — 1 in 
Ea. (|14|l . the initial [t ^ l/^z) isentropic drop in tem- 
perature can be written as 



(15) 



Here bp{t) satisfies the scaling equations for expanding 
hydrodynamic thermal clouds 8] in the presence of a 
mean field 0| , 



bl'hT 



5/3,4/3 



bV'b 



'b'^b'^' 



(16a) 



(16b) 



The Eas. (|16a|) and Ijl6b|l decouple for t <^ 1/wz since 
bz — 1. In this limit the radial scaling equation can be 
written as 

(^) -li^-m-i/b'At)]+^[i~i/bi{t)]. (17) 

We then substitute Ea. H17|) into H15|l and obtain to first 
order in [bp/ujpY the temperature reached at i = i*: 



2 / b„ 



Tn ^ 3 I Wo 



(18) 



t=t. 



We point out that in the limit of very elongated clouds 
Eq. (|17|l also represents the correct description for a fully 
hydrodynamic expansion. Then, we may write for the 
asymptotic expansion velocity in radial direction 



lim 



Kit) 



t^oo ujp Ip (0) V To 



Hence, comparing with the asymptotic value of Ea. H17f) 
we conclude that the following inequality should hold: 



1 < Tp/To < 3/2 + e 



(20) 



Returning to our experimental conditions wc empha- 
size that the duration of the hydrodynamic stage will 
be very brief because the instantaneous mean free path 
grows quadratically with bp in these elongated clouds. 



X{t)/Xo = bl{t), 



(21) 



as follows with Eqs.(^ and p4|l Js^. Roughly speaking 
is reached when the mean free path equals the radial 
size of the cloud. Therefore, a rough estimate for can 
be obtained by substituting A(i) — lp{t) into Ea. (|21() for 
t = tt,. With Eq.© this leads to 



bp{t^) ~ l/a)p7 



(22) 



As for t < l/u!p the radial size of a hydrodynamic cloud 
hardly differs from that of a coUisionless cloud. 



(23) 



bp{t) ^ {1 + ^yf\ 

we find with Eg. lf^ 

U ^ il/C^p)[{l/CupTaf - l]i/2 ^ 0.6 ms. (24) 

A self consistent estimate for our expansion model can 
be obtained by combining Eas. H23l) and H15|l for t — t^,, 



U ^ {l/Cjp)[{To/T,) 



3/2 _ -|ll/2 



(25) 



However, for this estimate the ratio T^/Tq should first be 
established experimentally. 
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B. Collisionless stage 

Once the expansion is ballistic {t > t^) the variance of 
the axial (i = z) and radial (i = p) velocity components 
of the expanding gas can be written as 

(vf) = (uf) + (wl) , (26) 

where Ui represents the thermal velocity components of 
the atoms and Wi the dynamic velocity components of 
the density distribution due to the expansion. 

At the start of the ballistic stage {t = t») the thermal 
velocity components can be associated with , 

m = keT, . (27) 

The dynamical velocities due to the overall expansion can 
be expressed as 

m (wf) = m (rf) = (5,/w,)'fcsTo = ^mll (28) 

Here we used the scaling property = {bi/bi)ri, with 
the Ti representing the position coordinates in the ex- 
panding cloud. Since for collisionless clouds the (vf) are 
conserved by the time the mean field has vanished, we 
may write 

m{vf)^kBT„ (29) 

where the Ti are effective axial and radial temperatures 
that may be associated with the asymptotic axial and 
radial expansion velocities si defined in Ea. (|ll|l . 

Substituting Eqs.l^B and ^ into Ea. (|^ we obtain 
for ^ < 1 

§^1^(1)1/ SIM' 

where the second term on the r.h.s. represents both the 
hydrodynamic and mean field contributions to the dy- 
namic motion at i = <C 1/tOz and the third term the 
mean field contribution to the dynamic motion for t > 
|39j. With Ea. (|17|l this results in the following relation 
between Tq, Tp and T, in expanding elongated thermal 
clouds 

^To + CTo^^T^+Tp. (31) 

This equation is valid for small mean fields provided <ti 
l/ujz and expresses the energy conservation during the 
expansion. It implies 

= n. (32) 
This also follows directly by writing in analogy to Eq. H3U|) 

taking into account that {bz/uiz)^ is negligibly small [40l |. 



V. THERMOMETRY 

The result (|32() shows that with our measurement of 
Tz we directly probe the temperature of elongated clouds 
at the end of the hydrodynamic stage. Knowledge of T^, 
allows to obtain with Ea. (|25|l a self-consistent result for 

within our expansion model. Using T^/Tq = 0.71(2) 
we calculate = 0.28 ms, somewhat smaller than the 
rough estimate (|24|l 

Rewriting H31|l we find an increase in the effective ra- 
dial temperature 

^ = |(l-i^) + e= 1.18(2). (34) 

Hence 15% of the increase in Tp is due to the mean field. 
Note that Ea. H34|) satisfies inequality (|20|l . Notice further 
that the value Tp = 1.37(6) obtained with Ea. l|M|) 
comes close to the value Tp = 1.35(6) /zK following di- 
rectly from the radial expansion. 

We found the fitting procedure for determining To, Tp 
and to be very sensitive for the detailed shape of the 
fit function. Choosing a simple gaussian reduces the esti- 
mated values for these temperatures by as much as 25%. 
However, this enormous systematic error does not affect 
the corresponding aspect ratios by more than a few parts 
in a thousand. We found more indicators that the aspect 
ratios are more accurately determined than the absolute 
values. Interestingly, we find for the aspect ratios stan- 
dard deviations of typically 1%, i.e. twice as small as for 
the absolute size 41] . This points to some form of error 
cancellation. Also the fit to the aspect ratio is somewhat 
better than those of the separate plots. 

Let us now turn to the results for the aspect ratios 
as presented in Fig.El Using Eqs.Q, (CDIl and (EH) the 
evolution of the aspect ratio can be expressed as 

I pit) ^ [(3/2 + 0- 1/2 (T,/ro)]i/W , . 

lz{t)~ [i + ^ + {n/To)u;^At~t*yV^^' 

where we presume t ^ 1/^^p as in Ea. ljlOII . By con- 
struction this form satisfies energy conservation. In this 
way our fitting function stays as close as possible to a 
fit to a solution of the scaling equations. Fitting Eq. H35() 
to the data using ^ = 0.03 and = 0.3 ms we obtain 
T^/Tq = 0.72(1). The fit is shown as the sohd hue in 
Fig.[21 The result agrees within experimental error with 
that obtained from the axial expansion data but the ac- 
curacy is slightly better. The method lacks the accuracy 
to extract ^ 42]. The dashed line in Fig.El corresponds 
to the collisionless limit of Ea. (|35|l : ^ = 0, = and 
n - To. 

Once Ea. H35|) is accepted, timc-of- flight information for 
a single expansion time sufRces for thermometry. The 
procedure goes in two steps. First we set ^ and t* equal to 
zero and use Ea. (|35|l to obtain a first estimate for T^/Tq. 
With Eq.^ Tp/To follows. After Tp is determined with 
Ea. (|10|) . we have an estimate for the absolute value To- 
Together with no, deduced with Ea. H12|l . this allows us 
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to calculate f and t^, . Iterating the procedure once yields 
all values within the limits of accuracy of the analysis. 
Choosing the expansion time sufficiently long {t ^ l/i^z) 
the results are very insensitive for the value of . 

Our estimates for the absolute values of Tq, Tp and Tj 
are sensitive for the detailed shape of the clouds. There- 
fore, deviations from the Bose shape will result in sys- 
tematic errors, in particular if the cloud shape changes 
during the expansion. Shape deviations can arise from 
the presence of the mean field. Also, inhomogeneous isen- 
tropic cooling as a result of the inhomogeneous density 
profile of our samples can give rise to deviations of the 
Bose shape. Further, it may be that our transforma- 
tion from optical density to column density gives rise to 
slight distortions of the cloud shape as a result of optical 
pumping or saturation of the detection transition. 

In our analysis we did not correct for deviations of the 
cloud shape from the Bose distribution. First of all be- 
cause under our conditions the mean field is weak and our 
fits of Eq.Q to the measured column densities look con- 
vincingly. Secondly, because shape deviations produce 
similar relative errors in all three temperatures. There- 
fore, they do not affect the conclusions and consistency 
of our analysis as long as the scaling approach remains 
valid. 



VI. CONCLUSIONS 

We studied the behavior of dense elongated clouds of 
*^Rb in the crossover from the collisionless to the hydro- 
dynamic regime. At our highest densities the mean free 
path is slightly smaller than the radial size of the cloud 
and the expansion is anisotropic. The expansion can be 
described by a two stage model in which the expansion 
is treated as purely hydrodynamic up to time t — 



and as purely collisionless beyond this point. We find 
that at the end of the hydrodynamic stage the tempera- 
ture has dropped substantially due to isentropic cooling, 
r*/To = 0.72(1). This reflects itself in an axial expan- 
sion that is substantially slower than expected for the 
collisionless case, = T*. In accordance with energy 
conservation the radial expansion is faster, Tp > Tq. The 
isentropic cooling is best determined from the aspect ra- 
tio. Although the mean field in the trap center is sub- 
stantial, L'nif(0)/fcTo = 0.23, it hardly affects the expan- 
sion behavior. Including the mean field in the analysis 
only affects the value obtained for Tg, with Tz and Tp by 
definition being unaffected. In our case the mean-field 
corrections are too small to be extracted with a fitting 
procedure but can be calculated accurately. It leads to 
systematic errors in Tq of order 3% if not included. 

Presently it is possible to study the case of strong inean 
fields by tuning to a Feschbach resonance |^ 0, IitL 
JJ., 43, 44J. It would be interesting to study the case 
where anisotropic expansions are to be expected, but the 
behavior of the system is dominated by the mean field 
rather than the coUisional hydrodynamics. 
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